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Abstract 

We investigate the effects of the atomic interaction in the Superradiant Rayleigh scattering from 
a Bose-Einstein condensate driven by a far-detuned laser beam. We show that for a homogeneous 
atomic sample the atomic interaction has only a dispersive effect, whereas in the inhomogeneous 
case it may increase the decay of the matter-wave grating. 
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I. INTRODUCTION 



The long coherence time of the Bose-Einstein condensate (BEC), now routinely produced 
in many laboratories offers the possibffity to study the eolleetive motion indnced by exter- 
nal radiation beams In particular, a single far-off resonance laser sent on an elongated 
BEC produces superradiant Rayleigh scattering 

Baa, 

generating coherent backscattered 
radiation and splitting the condensate into fractions moving at velocities differing by multi- 
ples of 2hk/m, where k = uo/c is the wave- vector of the laser incident along the symmetry 
axis of the condensate, u is the laser frequency and m is the atomic mass. Superradiant 
Rayleigh scattering from a BEC is the quantum analog of the collective atomic recoil laser 
(CARL) 5] in which the emitted radiation is not confined in a high-Q ring cavity [jj. The 
complete absence in a BEC of Doppler broadening due to thermal motion allows for a regime 
in which the atoms scatter a single laser photon and recoil with an extra momentum of 2hk 
in the direction of the incident photon. The number of scattered photons and the amplitude 
of the density grating resulting from the interference between the two atomic wavepackets 
with momentum difference 2hk are exponentially enhanced via the CARL instability j?| . In 
the absence of any mechanism of atomic dephasing, the process is sequential 0, 8|, with a 
complete transfer of atoms, after the superradiant process, from the original motional state 
with momentum p to a state with a momentum p + 2hk. However, in a real BEC several 
mechanisms contribute to the decay of the coherence between the two momentum states. 



Some of them are due to the decoherence induced by spontaneous emission p or phase 





diffusion [4]. The main characteristic of this kind of decoherence is irreversibility. Other 
mechanisms arise from inhomogeneous broadening, as those due to a finite size of the con- 
densate wavefunction, responsible for a broadening of the atomic momentum distribution, 
and from mean-field broadening due to the atomic interaction [10]. It has been recently 
suggested that the dephasing due to inhomogeneous broadening can be reversed applying, 
after the superradiant scattering process, a Bragg pulse of area ir inducing a superradiant 
echo and a further tranfer of the atoms to the final momentum state 3). In this paper we 
investigate the effects of the mean-field atomic interaction on the superradiant scattering 
process. 
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II. BASIC MODEL 



We consider an elongated Bose-Einstein condensate driven by a single laser incident along 
the positive direction of the symmetry axis z of the condensate. The laser is far-detuned 
from the atomic resonance, so that radiation pressure due to absorption and subsequent 
random incoherent, isotropic emission of a photon, can be neglected. In this regime, the 
atoms backscatter photons at frequency uj s and wave vector k s = uj s /c ~ k, recoiling with a 
momentum 2hk along the same direction of the incident laser beam. 

In a simplified ID description of the process along the axis z, the evolution of the matter- 
wave field ^(z, t) and of the dimensionless amplitude a(t) of the scattered radiation is 
determined by the following self-consistent equations: 

d 2 ^ 



a * e i(8St) _ c _ c _ 



^ + 2tt/3|^| 2 ^ (1) 



^ = gN f d0\*\ 2 e« e - st ') - m. (2) 

etc J 

where 9 = 2kz, a = (eoV/2hco s ) 1 ^ 2 E is the dimensionless electric field amplitude of the 
scattered beam with frequency u s , uo T = 2hk 2 /m is the two-photon recoil frequency, g = 
(Q/2A )(ujd 2 /2heoV) l l 2 is the coupling constant, VL is the Rabi frequency of the laser beam 
of frequency uj = ck, detuned from the atomic resonance frequency u by A = ui — u , d is 
the electric dipole moment of the atom along the laser polarization direction, V is the volume 
of the condensate containing iV atoms, 5 = u — u s and eo is the permittivity of the free 
space. The second term on the right hand side of Eq.(£Q) is the self-consistent optical lattice, 
resulting from the interference between the laser and the backscattered radiation, whose 
amplitude is amplified by the matter-wave grating described by the first term on the right 
hand side of Eq.(J2J). The matter-wave field \& is normalized to one, i.e. / d6\^(8, t)\ 2 = 1, and 
the last term on the right hand side of Eq. (0) describes the atomic interaction due to binary 
collisions, where (3 = AHka s N/mTi, a s is the scattering length and E is the condensate cross 
section. Eq.(j2J) has been written in the "mean-field" limit, which models the propagation 
effects of the light by replacing the nonuniform electric field by an average value and by 
adding to the equation a damping term with decay constant k ~ c/2L, where L is the 
condensate length and c is the speed of light in vacuum. 
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III. HOMOGENEOUS CASE 



If the condensate is much longer than the radiation wavelength and the density is uniform, 
then periodic boundary conditions can be assumed on 9 and the wavefunction can be written 
Fourier series 

^(9,t)=Y,c n (t)u n (9)e- mSt , (3) 

n 

where u n (9) = (1/ \/2ti) exp(in9) are momentum eigenfunctions with eigenvalues p z = 
n{2Tik). Using Eq.®, Eqs.(JTJ) and © reduce to an infinite set of ordinary differential 
equations, 

c n = -i5 n c n + #(a*c n _i - ac n+ i) - i(3^2c m cic* m+l _ n (4) 

m,l 

a = gNj2°nC* n+1 - na, (5) 

n 

where 5 n = n 2 u r — n5 and the dot indicates the time derivative. 



A. Two-level approximation 



Assuming that the only two momentum levels involved in the process are the initial level 
n and the final level n + 1, Eq.(jU) and (0) reduce to: 



c n+l 



W —I 

m —i 



S n + (3 (|c n | 2 + 2\c n+ i 
On+1 + /?(2| 



c n - gac n+ i 

c n+l + 9 a c r, 



a ~ gNc n c* n+1 — na. 

Defining 5* = c n c* +1 and W = \c n \ 2 — |c n+ i| 2 , we obtain from Eqs.tjHJ)-©: 

S = -i (A - /3W) S + gaW - 7 S 
W = -2g (aS* + c.c.) 
a = gNS — Ka, 



(6) 
(7) 
(8) 



(9) 
(10) 

(11) 



where A = 5 n — 5 n +i = 5 — uj r (2n + 1) and we have introduced a damping term in Eq.(|5|). 
to account for the decay of the coherence between the two motional states n and n + 1. We 
note that when the atomic interaction is neglected (/? = 0), A = is the Bragg condition 



of the scattering process, arising from momentum and energy conservation 



121. We observe 



from Eq. Q that the atomic interaction term has a dynamical dispersive effect on the Bragg 
resonance, proportional to the population difference W. In the linear regime, when a is 
still small and W « 1, the Bragg condition is A = (3, i.e. 5 = u r (2n + 1) + n a U/h, where 
n a U = n a A-Kh 2 a s /m is the chemical potential and n a = 2N/XT, is the atomic density. 

In the superradiant regime the field amplitude a can be adiabatically eliminated for times 
much longer than k~ 1 . In fact, let introduce the slowly varying variable Sit) = S(t)e %a<yt \ 
where a{t) At - f3 J * dt'W(t'), and let integrate Ea.flTTl): 

a(t) = a(0)e~ K * + gN C dt'S(t - £ ') e -ia(t-t')-«*'. ( 12 ) 

Jo 

If we assume that S and W do not change appreciably in a time k~ 1 during the superradiant 
process, i.e. if r sr ^> k~ 1 where r sr is a characteristic time for superradiance, then in Eq. (jl2|) 
S(t — t') w S(t) and a(t — t') m a(t) — [A — f3W(t)}t'. Performing the residual integration 
in t' and assuming t 3> k" 1 , we finally obtain: 

a(t) * 9NS(t) (13) 
k — i[A — j3W(t)} ' [L6) 

so that the field a follows istantaneously the atomic evolution. Combining Eqs.(jnj), (fTIHl and 

(fT3|) and defining I = \S\ 2 , we obtain 

f GW 1 

^ 2 { l + [ (A- W -^ ^ 

W = ^ (15) 

1 + [(A-/?W)/k] 2 1 ; 

where G = g 2 N/n is the superradiant gain. Defining the characteristic time as r sr = 1/G, 

it follows that the adiabatic approximation ()13|) is true for k 3> G, i.e. for k ^> g^/N. 

Furthermore, the two-level approximation is valid if G < u r . ^From Eqs. lfT^ and (fTB^) it is 

possible to derive the following analytical results: 

• In the linear regime, for W ~ 1, the threshold condition for superradiance is G > 
j{1 + [(A — /3)/k] 2 }, so that the only effect of the atomic interaction is a shift of the 
resonance from A = to A = (3. 

• Neglecting decoherence (7 = 0), Eqs. (|14|) and (|15|) admit the constant of motion 
AI + W 2 = 1. Writing 2\fl = sin0 and W = cos0, we obtain an equation for the 
Bloch angle 0: 

GSi °^ (16) 



1 + (A- ficoscpy /k 



Although Eqs.(fTl)]) can be solved exactly by quadrature, its solution can not be set 
in an explicit form. FigQ shows \a\ 2 /N, (a), and P n = \c n \ 2 = (W + l)/2, (b), as a 
function of u r t for different values of (3, obtained solving numerically Eqs.(HJ) and (j5J) 
for k = 20a> n gVN = 2u r and A = (3. The results are in excellent agreement with 
the numerical solution of the approximated Eq. pfij) . not reported in the figure. We 
note that the effect of the atomic interaction is only a broadening of the superradiant 
pulse, which still preserves the same area equal to ir, transfering completely the atoms 
from the initial momentum state n to the final momentum state n + 1. 

• It is easy to calculate the exact analytical solution of Eqs. (fi4")) and (fT3j) when (3 = 0. 
In fact, let introduce G' = G/[1 + (A/k) 2 ] and the new variables x = (W -W ) / '(1— W ) 
and y = 2y/l/(l - W ), where W = j/G' < 1. /.From Eqs. (HU) and (JT3J) it follows 
that x 2 + y 2 = 1. Introducing again the Bloch angle <fi defined such that x = cos 4> and 
y = sin 0, Eqs. (114)) and (fTHj) give the following equation for 0, 

= G"(l-^o)sin0,, (17) 

whose solution yields x(t) = — tanh[G'(l — Wo)(t — to)], where to = — ln[|*S'(0)|/ (1 — 
Wq)]/G'{1 — Wo) is the delay time. Coming back to the original variables we finally 
obtain: 

I(t) = 1^—^-) sech 2 [G'(l - W )(t - t D )} (18) 
W{t) = W - (1 - W ) tanh[G'(l - W )(t - t D )} (19) 

We note that the asymptotic value of the population difference is 2Wo — 1, so that the 
fraction of atoms left in the initial state after the superradiant process is P n = Wo = 
~f/G'. Measuring experimentally G' and P n it is possible to evaluate the decoherence 
rate 7 j^y]. 

IV. INHOMOGENEOUS CASE 

Let now consider the case in which the condensate is described initially by a wavepacket 
with a finite size oq = 2ka z and p z = 0. We have solved Eqs.Q and (0) for an initial 
Gaussian wavepacket of width ae = 25, k = I0u> r , g\/2N = cu r , A = and different values 
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of /3. The numerical integration of Eqs. and (J2J) is based on a finite-difference predictor- 



corrector scheme 



3- FigOshows the density distribution, p(6,t) = \^{6,t)\ 2 , for /3 = 



at different times, whereas figEl shows the corresponding momentum distribution p(p z ,t) = 
\^(p z , t)\ 2 , where \& is the Fourier transform of the wavefunction We observe that the 
superradiant process produces a condensate fraction moving with an average momentum 
p z = 2hk and a smaller momentum spread (see figE}. In the configuration space (see figHJ) 
we clearly observe the interference fringes when the two fractions overlap, whereas for longer 
times the recoiling atoms move away from the original condensate. 

FigJUshows \a\ 2 /N, (a), and P n , (b), as a function of ui r t for different values of (3, where P n 
is the population of the momentum state p z = n(2hk), calculated integrating the momentum 
distribution over an interval centered around p z = n(2hk) and of length 2hk. This can 
be done only if the momentum distribution remains narrower than the momentum level 
separation, i.e. if a Pz <C 2hk. 

We observe that contrary to the homogeneous case, increasing f3 the superradiant pro- 
cess becomes less efficient, decreasing the area of the superradiant pulse (see figE^a)) and 
increasing the fraction of atoms left in the initial momentum state (see figHJb)). This effect 
can be interpreted as due to a dephasing caused by a detuning from the resonance depend- 
ing on the atomic density. Each atom evolves with a different detuning from the resonance, 
resulting in a inhomogeneous broadening of the superradiant transition and a subsequent 
decay of the coherence between the atoms. Similarly to the photon echo jllj], it is expected 
that this dephasing may be partially reversed applying a suitable Bragg pulse or area n, at 
least for small values of (3. For larger or negative values of (3, it is expected that nonlinear 
effects are more important, and this regime will be the object of a future detailed analysis. 
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FIG. 1: Effects of the atomic interaction on the superradiant regime in the homogeneous case: 
\a\ 2 /N, (a), and population fraction P n of the initial momentum state, (b), vs. u r t, from the 
numerical integration of Eqs.Q) and © with k = 20uj r , g\^N = 2ui r , A = (3 and different values 
of f3. 
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FIG. 2: Evolution of the density distribution p(6,t) = \^(9,t)\ 2 vs. 9 = 2kz at different times, from 
the numerical integration of Eqs.Q and @ for an initial Gaussian wavepacket of width ag = 25, 
k = 10u r , gV2N = u r , A = 0, and /3 = 0. 
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FIG. 3: Evolution of the momentum distribution p{jp z ,t) = \^{p z ,t)\ 2 vs. p z (in units of 2%k) at 
different times and for the same case shown in figEl 
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FIG. 4: Effects of the atomic interaction on the superradiant regime in the inhomogeneous case: 
\a\ 2 /N, (a), and population fraction P n of the initial momentum state, (b), vs. uj r t, for the same 
initial conditions and parameters of the case shown in figEl and different values of (3. 
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